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Abstract 

We propose a method for a QCD based calculation of one-particle 
inclusive decays of the form B —> DX or B ^ D*X. It is based on 
the heavy mass limit and a short distance expansion of the amplitudes, 
which yield a power series in the parameter 1/M^ for the spectra and 
in AgcDirib/imb — rricY the rates. We study the leading term of 
this expansion for the case of the semi-leptonic decays B —> DXi'^v. 



1 Introduction 


Over the last ten years significant progress has been made in the theoretical 
description of heavy flavour decays [|l|. The application of the l/mg ex¬ 
pansion [rriQ being the mass of the heavy quark) allows us to perform QCD 
based calculations, which in some cases yield model independent results. The 
additional symmetries of the infinite mass limit, the so called heavy quark 
symmetries P], reduce the uncertainties due to unknown hadronic matrix 
elements significantly, and corrections to this infinite mass limit have been 
studied extensively using the framework of Heavy Quark Effective Theory 
(HQET) §. 

The heavy mass expansion has been applied to various classes of decays. 
As far as exclusive decays are concerned the main progress has been achieved 
for semi-leptonic decays, while exclusive non-leptonic decays still have not 
simplified through the heavy mass limit. 

The other side are the fully inclusive decays, the rates of which can be 
obtained as a power series in l/mq by means of an operator product expan¬ 
sion (OPE) and subsequent application of HQET IQ. Here semi-leptonic as 
well as the non-leptonic processes may be described, allowing us to compute 
lifetimes and branching ratios. The pattern is well reproduced by the l/mq 
expansion, although some open problems remain 

Up to now no attempt has been made to apply similar methods to one- 
particle inclusive decays, such as B ^ or B ^ DX and B D*X. 

Obviously the standard method as in the inclusive case does not work in a 
naive way, since in the final state a D or a D* is projected out, and the same 
set-up as in the inclusive case will not work. 

In the present paper we propose a method which allows us to compute 
one-particle inclusive rates, based on QCD. The main ingredients are similar 
as in the fully inclusive case. In the next section we shall describe the method 
and then discuss its application to semi-leptonic decays. 

2 Description of the Method 

We shall consider first decays of the form B DX (i.e. a 6 —>■ c-transition) 
and thus study the expression 

G(M=) = Y,\{B(pB)\H,„\D(p[,)X}f i2w)VipB -PD-Px) ( 1 ) 
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where |X) are momentum eigenstates with momentum px and H^ff the 
relevant part of the weak Hamiltonian. The function G depends on the 
invariant mass = {ps —PoY of fh® state |X) which ranges between 


0 < < (niB — rnoY 5 


( 2 ) 


where we have neglected the pion mass as well as the lepton masses. This 
function G is related to the decay rate under consideration by 

dT{B ^ DX) =G{M^) (3) 

2uib 

where is the phase space element of the hnal state D meson. 

The region close to 0 is dominated by a few resonances (the vr 

and p states in the non-leptonic case), and away from this region one can 
expect duality to hold. In particular, this should be true in the limit in which 
TTib, rUc —»• 00, since in almost all available phase space we have 3 > ^qcd- 
In technical terms this means that we are going to set up a short distance 
expansion for the quantity G{M^). The procedure is similar as the one for 
inclusive decays, we write 

X 

( 4 ) 


and make use of the fact that rric and mf, are both large scales. We make 
these scales explicit by redehning the heavy quark helds in H^ff by 

b{x) = c(x) = (5) 

where the velocities are dehned as ps = and po = m^v'. 

Inserting this yields 

G{M‘^) =J2[ ( 6 ) 

{B(v)\Heff(x)\D(v')X){D(v')X\H,„{0)\B(v)), 

where H^ff is obtained from H^ff by the replacements b ^ by and c —>■ Cy^. 
Equation shows that the large momentum entering the game is rUbV—rricv'. 
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The next step is a short distance expansion of the matrix element appear¬ 
ing in yielding a power series in inverse powers of the large momentum 
M = nibV — rricv' 


G{M^) = Y,E^‘"'MWps)\ot>\B{ps))U. ( 7 ) 

n=0 i 


The operators depend on the final state Z)-meson and are the analogue 
of the production operators as they appear in heavy quarkonia production 
IP or in one-particle inclusive production in annihilation p. They are 
local and have the generic structure 

(^ 1 "^ = J2\-^.'rK]\D{v')X){D{v’)X\[Krc,^ , ( 8 ) 

x 

where T^'^ denotes a combination of Dirac matrices and covariant derivatives. 

The matrix elements of the between static B meson states are univer¬ 
sal functions of the velocity product v ■ v'. We shall not give a detailed proof 
of factorization of the matrix elements into long and short distance contribu¬ 
tions, rather we are aiming at a phenomenological analysis of the one-particle 
inclusive semi-leptonic decays. We remark that the method is not as rigorous 
as in the case of fully inclusive decays, where the heavy mass expansion is 
derived by an operator product expansion. However, fig.|5 makes the argu¬ 
ment plausible. For large enough the large momentum flows through the 
state |df) and we assume that parton-hadron duality holds for this part of 
the diagram, and hence we can compute this part in perturbative QCD. 

The dimension of in (|^) is n-|-6 and hence is of the order 

M. The leading term of the expansion involves dimension 6-operators and 
we shall discuss in the present paper only this contribution. If we consider 
only this leading term, we may even replace the operators and Cc' by static 
HQET quarks. In the following this replacement is understood. 

The corrections to be expected can easily be estimated. Since the full 
momentum transfer is Q = M + k where k is the sum of the residual momenta 
of the heavy b and c quark, the corrections to the leading term originate 
typically from 

Q^ = M^ + 2M-k + 0 {K%cd) =M^(^ + (9) 
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Figure 1: Illustration of the short distance expansion of 


and hence the corrections involve typically matrix elements of the form 

5^(S(n)|[c,.r6J|D(n')X)(D(n')X|(M-^D)[6,rtc,,]|5(n)) (10) 


where D is the covariant derivative of QCD. This matrix element will be of 
order M ■ vAqcd and consequently the method works as long as 


2M • n, , , 

—j^Aqcd < 1 • ( 11 ) 

and M-v are not independent variables, since they both can be expressed 
in terms of the velocity product v ■ v'. Eliminating M ■ v one obtains 


Aqcd 

mb 


2 2 
AP 



< 1 , 


( 12 ) 


and hence it is obvious that the expansion breaks down for very small 
Here again a similar situation occurs as in the inclusive semi-leptonic decays, 
where the endpoint region may be described in terms of a shape function. 

On the other hand one may ask wether the short distance expansion works 
at all, and thus it is instructive to insert the maximal value for which is 
possible in a decay. One hnds that the parameter 

‘^Aqcd 
mb - me 
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should be small compared to unity. Inserting the pessimistic value Aqcd = 
500 MeV one hnds that this parameter is about 1/3, which justihes our 
approach for the spectra at least close to maximal M^. In order to get the 
total rates an integration over the phase space of the D meson has to be 
performed. The details depend on the process under consideration, but the 
typical size of the corrections can be estimated by computing some arbitrary 
phase space average. Choosing a pase space measure as 

= 2mh dp — , 

\/[{mc + rabY — M^] [{rric — rribY ~ 



which yields very simple integrals we hnd 

M -Vx _ J ^ mbAqcD ^ 1 

' M2 ^ ~ {mb - mcY ^ 4 


(14) 


justifying the short distance expansion also for the rates. 

The second type of corrections are the QCD radiative corrections which 
can be computed systematically. They will be of the order as{M‘^) and hence 
will be small enough to be treated perturbatively. As usual, the logarithms 
as{M‘^) In(M^) can be resummed by renormalization group methods; for the 
leading terms this will be done in section |^. 

In the present paper we study only the leading term of the expansion and 
focus on applications to weak interactions. In this case we need to consider a 
matrix element of a dimension-six operator involving the left handed currents. 

We consider the leptonic case in some detail; inserting the well known 
effective Hamiltonian for semi-leptonic decays we hnd 

G(M2) = (15) 

5^(P(n) I [c,,7'^(l - 75)M l^(^'')^) - 75)c.^] |i?(^)), 


where is a tensor originating from contracting the lepton helds in the 
effective Hamiltonian. This tensor only depends on the vector M and hence 
has the form 


P^,(M) = A{M^){M^g^, - M^M,) + B{M^)M^M, (16) 
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Neglecting the lepton masses, we obtain at tree level 

A{M^) = ©(M^) and = 0 . (17) 

Stt 

Using 

- 75)c^/ = {mb - mc)Kc^> - {mb + mc)K'^bCv' (18) 

we can write the leading order contribution as 

G(M^) = ^\Vcb\'^AmBmD [{rriB - mDfTls{v ■ v') (19) 

DTT 

+ {mB + mo^Vpi'^ • "^0 “ M‘^{r]v{v ■ v') + t]a{v ■ U ))] 
where we have dehned non-perturbative matrix elements as 
AmBmDris{v ■ v') = '^{B{v)\[c^iK]\D{v')X){D{v')X\%c^,]\B{v)) (20) 

-AmBmDr]p{v ■ v') = '^{B{v) \ [c^''ybK]\D{v')X) {D{v')X\[b^'y5C^>]\B{v)) 

X 

AmBmDriv{v ■ v') = '^{B{v)\[c^i')%^]\D{v')X){D{v')X\%'^^Cy>]\B{v)) 

X 

AmBmDr]A{v ■ v') = ^( 5 (n)| %i'y^^'-f^h^]\D{v')X){D{v')X\%'y^'-f^c^,]\B{v)). 

X 

Once radiative corrections are taken into account, these operators mix with 
the corresponding operators where the b and the c quark are coupled to a 
color octett: 

AmBmDPs{v ■ v') = \ [c,n]\D{v’)X) {D{v’)X\[KT^c,,]\B{v)) ( 21 ) 

X 

-AmBmDPp{v ■ v') = ^(B{v)| [c„/75T“6„] |Z)(U)X) (i 7 (U)X| [6^75T“c„/] | 5 (t;)) 

X 

AmBmnPv{p ■ v') = ^(-B{v) | [c„/7'"T“6„] |Z1(U)X) (Z1(U)X| [6^7^T“c„/] |5(t;)) 

X 

AmBmDPA{v ■ v) = 

J2{B{v) I [c,'Y^,T%,]\D{v')X) {D{v')X\[K^^^,T^c,,]\B{v)) 

X 

Note that we are using parton-hadron duality and thus |X) is expressed 
in terms of QCD degrees of freedom and thus the matrix elements p* are 
nonvanishing. 
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3 Renormalization Group Improvement 


Under renormalization the matrix elements become scale dependend 

quantities rjiiy ■ pi{v ■ v',p). We chose to construct the short distance 

expansion at an intermediate scale rh = l/2{mb + mc), where both the b- and 
c-Quark are treated as static helds described by the HQET. However, the 
typical scale of the hadronic matrix elements is a low hadronic scale p = A. 
Using the renormalization group the matrix elements can be scaled from the 
matching scale rh down to the low scale A. In our numerical analysis we use 
^^(A) = 1 , where as{p) is the one loop expression for the running coupling 
constant. 

To this end one has to renormalize the operators [] 


X 

= Y}^.'Tn]\D{v')X){D{v')X\[KV{r^c,,] (23) 

JV 


where Tj = 1 , 75 , 7 ^, 7 ^ 75 . Because of heavy quark spin symmetry mixing 
occures only between singlett and octett operators corresponding to one spe- 
cihc Dirac structure Tj. That means a basis closing under renormalization is 
given by and (9® for every individual i. 

The mixing properties of the operators translate into that of their matrix 
elements. Therefore we can formulate the renormalization group equation 
directly in terms of the r]i{y ■ v\ p) and piiy ■ v\ p) as follows: 

■ v\ p) = 7ii?7i(n ■ v\ p) + iispi{v ■ v\ p) 

-j^pi{v ■ v', p) = 'y%iPi{v ■ U, p) + 788P^(^^ • v\ p) 
amp 

Since we restrict ourselves to the leading logarithmic approximation, it suf- 
hces to know the one loop anomalous dimensions. These are given by the 

’^Note that at least to one loop order the renormalization properties of these operators 
are identical to those of the operators 

Df) = [c,,TX][KT,c,,] = %,T,T-K][KT,T-c ^,]. (22) 

since the UV behaviour is independent of the states, including the D appearing in the 
final state. 


(24) 
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Figure 2: Feynman diagrams contributing to the one loop anomalous di¬ 
mensions. The blob represents generically the operators corresponding to 
the rji and p*. 


divergent parts of the Feynman diagrams shown in hgure ^ supplemented by 
wave function renormalization of the heavy quark helds 


c\ \ rv 

711 = ■ v') 718 = -{^)2K{v ■ v') 


IT 


TT 


781 = (—)^(^ - ■ v') 788 = {—)^K{v ■ v') 

TT 1 Ni TT Nr 


(26) 


where Nr = 3 is the number of colors and 


K{v ■ v') = 1 — V ■ u'Re[r(u ■ u')] r{z) = 


\n{z + — 1) 


(26) 


The function r{v ■ v') typically appears in the anomalous dimensions of ve¬ 
locity chanching heavy quark currents 0. 

In our case only the real part of r(u • v') shows up in the anomalous 
dimensions, since the corresponding Feynman amplitudes contribute to the 
forward matrix element of a hermitian operator which has to be real. Note 
that individual Feynman diagrams develop imaginary parts which drop out 
in the sum. 
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Solving we express the matrix elements T]i, pi at the scale m in terms 
of their value at an arbitrary scale p: 


Ei{v ■ v') = 7]i{v ■ v\m) = Cii{v ■ v\p)r]i{v ■ v',p) + C'i8(n • v\p)pi{v ■ v\ p) 

Ri{v ■ v) = Pi{v- v', m) = Cfii{v ■ v, p)piiv ■ v', p) + ■ v', p)pi{v ■ v, p) 

(27) 


The coefficient functions Cij{v ■ v',p) are given by 

Cn{v-v',p) = + {1 - ^)C{vv', p) 

Cisiv ■v',p) = v', p)^ 

Csi{v ■ v', /^) = ^(]^ - 1) /^) - 

C88(i^ • n', /i) = 1 + ^ • v', p) - 


(28) 


where 


C{v-v\p) 


Nc 

/ a,(p) \ 2/3o 
[as(rh)J 


K(v-v') 


with /3o = (33 — 2Nf)/12 and Nf = 3 for three active quark flavours. 

Note that in the case of semi-leptonic decays only the functions Ei{v ■ v') 
are needed, since in LLA there are no octett contributions at the matching 
scale. 


4 One-particle Inclusive Semi—leptonic Decays 

We shall hrst try to understand the data on the decays B —>■ EXt^v. In 
order to do this we need to have some idea about the matrix elements pi 
and Pi {i = S', P, V, A) which are dehned in (^) and (PT|). We shall work to 
leading order in the 1/M expansion and hence identify rric = rriD = mD* and 
nih = niB. 

We are aiming at the energy spectrum of the D meson in the one-particle 
inclusive decays of the type B DXi~^u. The rate is obtained by integrating 
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over the phase space of the D. Taking into account renormalization one gets 


dV 

dy 


2m B 

3 _ 

127r3 




Vcb\^m]^^/y'^ - 1 [(ms - mofEsiy) 


+ {mB + mDfEp{y) - M‘^{Ev{y) + EA{y))\ 


(29) 


where y = v ■ v' and the Ei {i = S, P, V, A) are the renormalization group 
invariant combinations of the rji and pi 


Eiiy ■ v') = Ciiiy ■ v',p)T]i{v ■ v',p) + Cis{v ■ v',p)pi{v ■ v',p). (30) 

The Wilson coefficients Cn and Cig have been given in (P5|). 

To get some expression for the functions p, and pt {i = S, P, V, A) we 
hrst observe that ai v -v' = 1 the inclusive rate is saturated by the exclusive 
decays into the lowest-lying spin symmetry doublet D and D*. Furthermore, 
at this point only the pi contribute, since Ci%{v ■ n' = 1) = 0. The D* 
subsequently decays into D mesons and thus at v ■ v' = 1 the sum of the 
exclusive rates for B —> and B —^ is equal to the one-particle 

inclusive semi-leptonic rate B —> D£~^i>X, which is again equal to the fully 
inclusive rate B In other words, at this point there are no decays 

into other charmed hadrons than D mesons. 

Off this point things become more complicated. However, as far as the 
total rates are concerned, still the exclusive decays B Dt^v and B —> 
saturate the fully inclusive rate B Xc£^p at a level of about 70%. 
Since the D* decay all into D mesons, it is certainly a good starting point to 
approximate the p* by something one obtains from the sum of the exclusive 
decays. In other words, we shall express the pi in terms of the Isgur-Wise 
function |^. 

The approximation we are going to use corresponds to some kind of fac¬ 
torization assumption formulated for G{M‘^). The functions pi are dehned 
by the matrix elements (^) and we shall approximate these matrix elements. 
However, as with the usual factorization, our approximation is not a scale 
invariant concept, and hence we have to dehne, at which scale it should hold. 


10 







( 31 ) 


At a small hadronic scale A we replace in (POP 

I [c,>rA]\D{v')X) {D{v')X\[KriC,,]\B{v))l=A 

X 

{B{v)\c^Xiby\D{v'))\f,=A{D{v)\KTiCy'\B{v))\^=A 
+ (5(r;)|c,ir,6,|5(r;0>^(5*))U=A(5(r;0r(5*)|6,r,c,i|S(r;))|^=A 

Y{D*) 

where Y{D*) is dehned by D* —>■ DY{D*), i.e. Y{D*) is either a pion or 
a photon originating from a D* decay. In the following we shall call this 
replacement factorization, since it is closely related to the factorization as¬ 
sumption known from non-leptonic decays. We get, again schematically 

I [cyYA]\D{v')X) {D{v')X\[byr,Cy,]\B{v))\,=A (32) 

{B{v)\cyYiby\D{v'))\^=A{D{v')\byTiCy>\B{v))\^=A 
+ Y,(B{v)\cyYA\D*{v', e))\,=A{D*{v',e)\byr,Cy,\B{v))l=A 

Pol 

DY{D*)) 


where the sum runs over the polarization states of the D*. In (|3^) we have 
used the narrow width approximation for the D* in the intermediate state. 

The matrix elements appearing in the factorized expression (^) can all 
be expressed in terms of the Isgur-Wise function: 


{B{v)\cyfTiby\D{v'))\^ = ^^msmDTr{ 75 (l ^)ri(l -h /)75}^(^^ • ii) 

{B{v)\cyYiby\D*{v\ e))|/, = ^i/mim^Tr{75(l ■ n',p) 

(33) 


From this we get 


T]i{v ■ v\^i) 


\X{vY)\^ 


Ci{v ■ v') + c*(n • v')Bt{D* DX{D*)) 


where 

Ci{v ■ n') = ^|Tr{75(l + ^)ri(l + /)75}P 
Id 

Ci{v-v') = ^^|Tr{75(l + |^)ri(l + V")/}r 

Pol 


(34) 
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and X [v ■ v') is the renormalization group invariant combination 


A '(v ■ v') = C:i(v ■ v', ■ v', li). 


( 35 ) 


The Wilson coefficient renormalizing the > c„/ current is 

known to two loops, but since we computed Cn and Cig only to one loop, it 
is sufficient to insert the one loop result 


Cs{v ■ v',n) 


as{m)) 


(36) 


where 


7hhiv ■ v') 




If 

V r[v • V ^ 


(37) 


The factorization assumption yields expressions for the matrix elements 
rji at the small scale A, but it does not tell us anything about the color octett 
contributions p*. It is well known that factorization should hold in the limit 
Nc ^ oo. This fact is indeed reflected in the W^dependence of the Wilson 
coefficients, since 


lim Cis = lim Cgi = 0 

Nc^OO Nc^OO 

lim Cgg = 1 (38) 

Nc^OO 

lim Cn = ICsp 

and thus the dimension 6-operators renormalize as products of dimension 3- 
currents and factorization becomes scale independent. This does still not tell 
us much about the pi, but a natural assumption is that they are of the order 
l/Nc and hence we shall take pi to be constant with pi{v-v', p) = 1/Nc- This 
simple ansatz, ignoring a possible dependence on v ■ v', does not introduce 
large uncertainties for the one-particle inclusive semi-leptonic decays, since 
the Pi are only induced through radiative corrections. 

In the following we shall consider the decays of the and the 5°, both 
of which contain a b quark undergoing a semi-leptonic transition b —> cT+z/. 
In the heavy mass limit for the c quark the hnal states involving a c quark 
{D^ or states) are suppressed, since this would involve a cc pair creation. 
To leading order in 1/mc the possible decays are thus 

B+ ifi+uX B+ D-tvX (39) 

^0 ^ ^ D-i+uX . 
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Since many of the D mesons originate from D* decays we have to take into 
account the relevant branching ratios of the D* mesons into the D mesons 
of different charge. We use o 

Br(D*- ^ Tfx) ^ 68 % Br(D*- ^ D-X) ^ 32 % ( 40 ) 

D^X) ^ 100 % 


and hence we can have the following decay chains 


B+ 
5 ° - 

- 



Br(D*o- 

-^D°X) 

5+ - 

D^£+uX 

D*-£+iy 

Br(L>*-- 

^D-X) 

5° - 

D-£+vX 

D*-£+iy 

Br(D*-- 

-*D°X) 

B° - 

^ D^+uX , 


(41) 


where the arrow indicates that - in addition to the direct decay channel 
B —> D£~^u - the exclusive mode on the l.h.s. contributes to the one-particle 
inclusive rate on the r.h.s. weighted with the branching ratios (^). 

We shall label the Ei for the different decay modes (|5^ with a superscript 
indicating the initial B and the hnal D meson. Taking into account the D* 
branching ratios (p][ ) we arrive at the following expressions for the Ei involved 
in the B~^ decays: 


Er'^iv) 


+ 2)(v + 1)|A'(J/)P + C,8(»,A)^ (42) 


In the B^ decays we have to take into account the D* branching ratios as 
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1 


ET^{y) 

Ef^iy) 

Ef^iy) 


Ef^'iy) 

Ef^'iy) 

E^-iy) 


E 


B°D- 

A 


(v) 


Ci,(9.A)i 

i\c 


^ + C'lsfe. A) 

^^Br(Z)- ^ S"A')i(!/^ - l)|A'(s)P + C,,(y, A) 
^ B°A)i(s + 2)(s + l)|A(s)|^ 


Nc 

1 


+Cl8(l/, 

I\c 


(43) 


|f;^ife + lf|AMP + C.,fe.A)i (44) 

£>-A)i(j,^ - 1 )|A(s)P + C,8(j,,A)^ 
(i(, + 1) + Br(D- -. D-X)l(y^ - 1)) |A(,)P 

+Ci8(!/.A)^ 

i\c 

~ ci(yA) ^r(0’- - D-X)^(y + 2)(y+ 1)|A(»)|2 
EEis(y, E)-—. 


Note that at v ■ v' = 1 we have simply the sum of the exclusive channels 
B —> Di^u and B —>• D*£~^u, where the D* component is weighted with the 
appropriate D* branching ratios, since here Cu = C 3 = 1 and Ci^ = 0. Off 
the point v -v' = 1 we still have C'ii/(C' 3 )^ ~ 1 but there is also an additional 
contribution from the octett contributions pi. As we shall see, this additional 
contributions are consistent with the data, despite of our crude estimate. 

Inserting these lenghty expressions into the master formula (P^ one 
obtains expressions for the one-particle inclusive energy spectra of the D 
mesons. To leading order in l/rric the energy of the D* meson is equal to the 
energy of the D meson originating from the decay D* —^ DX since X is soft 
of the order l/rric- 
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Figure 3: Decay spectra of the one-particle inclusive decays. Solid line: B — 
D~i^iyX, dotted line: B dashed line: B —>■ {D~ + D^)i^i'X, 

dashed-dotted line: 5° — {D~ -|- D*~)i^p. 


In order to actually obtain numbers one needs the Isgur Wise function 
as an input. A good £t to the experimental data is obtained already with a 
linear function, which is fitted to the renormalization group invariant X [y) 


xiy) 




with 


a = 0.84 
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( 45 ) 


In hgure ^ we plot the spectra of the D meson for the combined rates 








Mode 

Br (theory) 

Br (data from 11 ) 

B 

D-i+uX 

2.3% 

(2.7 ±0.8)% 

B 

'^i+uX 

6.9% 

(7.0 ± 1.4)% 

5° - 

^ D-i+u 


(1.5 ±0.5)% 

5° - 

^ D-*i+u 


(4.68 ± 0.25)% 

B 

D**iu 

(3.5 ±0.6)% 

(2.7 ±0.7)% 

B 

non-D£~^i> 

(0.6 ± 0.4)% 



Table 1: Comparison of our results with data. To get branching ratios, we 
used tb+ = tbo = 1.55 ps. Here D** denotes any hnal state with a D meson 
which does not come from the exclusive decays listed in row three and four. 
The last two rows are commented in the text. 


and compare it to the sum of the exclusive decays ^ D and —>■ 
D*-i+u . 

One may also integrate the spectra to obtain a total rate for the one- 
particle inclusive semi-leptonic processes. In table |l| we compare the rates 
we obtain from our approach with the experimental data from 


Table |I] and also hgure exhibit a few interesting features. First of all the 
experimental data are well reproduced. Furthermore, although we have used 
the assumption (|3^) our result is not simply the sum of the inclusive decays 
B —>• and B since (^2]) is a scale dependent statement. We 

assume that (^) holds at the small scale A; running up to the matching scale 
m yields a signihcant contribution from gluon exchanges. We interpret these 
contributions as B ^ where D** now stands for all D-meson hnal 

states, which do not originate from B —>■ or H —>■ Although 

the ansatz for the octett matrix elements pi is extremely simple, we obtain 
a reasonable number, namely Br(i? —>• Ri 32% x Br(B Xci~^u) 

where we use Br{B Xci'^v) = (10.4 ± 0.4)% from |^. From hgure ^ it is 

obvious that the Z)**-contribution vanishes at u • u' = 1 as required by the 
heavy quark limit. 

The last row of table ^ gives the branching ratio for decays which do not 
have a D meson in the hnal state, rather some other charmed hadron. The 
only other ground state hadron is a A^ so this should be the branching ratio 
for B —> for which we obtain Br(H —> = 6% x Br(H —> 
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Xg£+z/). This is what one would expect on the basis of the naive reasoning 
that a heavy quark hadronizes into a baryon with a branching ratio of about 
ten percent. 


5 Conclusions 

Exclusive semi-leptonic as well as fully inclusive decays of heavy hadrons 
have a well established basis in QCD. While in the former case it is the 
heavy mass limit of QCD, formulated as an effective theory (HQET), in the 
latter case it is the heavy mass limit combined with parton-hadron duality, 
formulated as an operator-product expansion. 

On the other side there are the exclusive non-leptonic decays, where no 
theoretically solid basis for a calculation of branching ratios exists. However, 
these decays are of prime interest with respect to CP violation and the de¬ 
termination of the CKM matrix. In this held the heavy mass limit has not 
brought any signihcant progress. 

In this work we have set up a QCD based description for one-particle 
inclusive decays. The basic ingredients are the heavy mass limit and a short 
distance expansion. We obtain operators similar to the ones describing heavy 
quarkonia production or one particle inclusive processes. We have formulated 
this method for decays of the type B DX where X in principle can be 
any state. 

We have applied this method to one-particle inclusive semi-leptonic de¬ 
cays of the form B DXi~^u, studying the leading order in the operator- 
product expansion. Higher order terms are suppressed by inverse powers of 
a large scale related to the heavy quark masses. To leading order, all these 
decays are parametrized in terms of eight functions rji and pi which depend 
on the velocities of the B and the D meson. 

The main problem is to obtain these non-perturbative functions rji and pi 
and we employed the fact that the inclusive semi-leptonic decays are dom¬ 
inated by the two channels B Di~^u and B D*i^u. Using this as a 
starting point we may obtain four of the unknown functions (the pi) in terms 
of the Isgur-Wise function with a well motivated factorization-ansatz. The 
remaining four (the pQ are suppressed by powers of as as well as by factors 
1/Nc. 

Estimating the four functions pi to be of the order I/Nq we are able 
to describe the features of the one-particle inclusive semi-leptonic decays. 
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In particular, QCD radiative corrections induce a relatively large amount 
of decays which originate not from the exclusive modes B —^ and 

B —> D*£^u. This is in accordance with the experimental data giving us 
some conhdence in our method. 

The approach suggested in the present paper opens the door to a QCD 
based description of one-particle inclusive processes; it is not limited to semi- 
leptonic decays. In particular, the non-perturbative functions rji and pi are 
universal and should also describe other one-particle inclusive processes. 
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